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Abstract 

Background: Sulphadoxine and pyrimethamine are anti-folate drugs tliat sliow synergistic anti-malarial effect. 
Point mutations in dihydrofolate reductase {dhfr) and dihydropteorate synthatase [dhps) cause anti-folate drug 
resistance phenotype in human malaria parasites. This study presents pattern of point mutations in dhfr/dhps genes 
among Indian sub-continent. 

Methods: Microscopically diagnosed one hundred Plasmodium vivax field isolates were collected from five widely 
separated geographical regions of India. Dlifr and dlips genes were PGR amplified and sequenced. Previously 
published mutations data were collected and analyzed using Chi square test to identify geographical cluster of 
mutant/wild type genotypes. 

Results: Sequence analysis revealed single (S58R), double (S58R/S117N) and quadruple (F57L/S58R/T61M/S1 17T/) 
point mutations at dlifr and single (A383G) to double (A383G/A553G) mutations at dlips in P. vivax field isolates. 
In addition, three new mutations were also observed at dlifr. Both, dhfr and dhps genes revealed tandem repeat 
variations in field isolates. Dhps revealed very low mutation frequency (14.0%) compared to dhfr (50.70%). 
Comparative analysis revealed a progressive increase in frequency of quadruple mutant dhfr genotype (p < 0.001) 
within five years in north-eastern state (Kamrup, Assam). Frequency of dhfr genotypes revealed three distinct 
geographical clusters of wild (northern India), double mutant (southern India), and quadruple mutant (north- 
eastern and island regions of India) on the Indian sub-continent. 

Conclusion: Study suggests that SP may be susceptible to P. vivax in India, except Andaman and north-eastern 
state. The distinction of geographical regions with sensitive and resistant parasite phenotypes would be highly 
useful for designing and administering national anti-malarial drug policy. 



Background 

Malaria is a life-threatening ancient parasitic disease and 
causes 250-500 million clinical episodes and nearly one 
million deaths annually [1]. Among the five human 
malaria species, Plasmodium falciparum is the most 
severe form, causing malignant malaria globally, while 
Plasmodium vivax is the most widespread species out- 
side Africa, causing huge morbidity and can be severe 
and fatal [2-7]. 

The worldwide spread of chloroquine (CQ) resistant 
strains of P. falciparum has led to use of sulphadoxine- 
pyrimethamine (SP) as the first-line anti-malarial drug 
in Southeast Asian countries. Sulphadoxine and 
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pyrimethamine sequentially inhibits dihydropteroate 
synthase (DHPS) and dihydrofolate reductase (DHFR) 
enzymes respectively in the folate biosynthesis pathway 
resulting synergistic anti-malarial effect [8]. Parasite has 
overcome the effect of SP by evolving point mutations 
in the respective genes encoding enzymes involved in 
the folate biosynthesis pathway. The mutated DHPS and 
DHFR enzymes have reduced binding affinity with SP 
drug and thus parasite survives in the presence of drug 
[9,10]. 

In India, resistance to CQ was reported for the first 
time in 1973 in P, falciparum from north-eastern states 
[11] and later it spread throughout the country [12]. To 
overcome CQ resistance problem, in 1982, SP was 
employed as first-line anti-malarial therapy for treat- 
ment of falciparum malaria in areas with >25% CQ 
resistant level, complicated malaria case, and higher 
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malaria endemicity [12]. Currently, according to the 
national anti-malarial drug policy, artimisinin-based 
combination therapy (ACT) (Artesunate+SP) is being 
used in most of the malaria endemic regions [12]. 

Reduced susceptibility to CQ in P. vivax was for the 
first time observed in 1989 in Australian soldiers 
returned from Papua New Guinea [13]. Later, several 
cases from Papua New Guinea, Indonesia, New Guinea, 
Brazil and India were documented [14-19]. Recent stu- 
dies from P. vivax ex-vivo maturation experiment 
showed reduction in the susceptibility to CQ in South- 
east Asian counties [20,21]. This information indicates 
gradual increase in CQ resistance cases of P, vivax, 

India contributes more than 78% of total malaria cases 
of Southeast Asia and P. vivax accounts for more than 
50% of annual malaria cases [22]. Point mutations in 
P, vivax dhfr had been documented from different parts 
of India [23-25], however mutation data of dhps is yet 
to identified in order to understand the molecular 
epidemiology of anti-folate resistance. Therefore, identi- 
fying information about presence of anti-folate drug 
resistance related point mutations {dhfr/dhps) in 
P, vivax from Indian sub continent would be highly 
helpful to understand the global pattern of anti-folate 
drug resistance. This study aims to identify point muta- 
tions in dhfr/dhps and the spatio temporal pattern of 
anti-folate drug resistance in Indian sub-continent. 

Methods 

Study sites and sample collection 

Blood samples were collected from five widely separated 
geographical regions of the Indian subcontinent namely 
Delhi (2005); Chennai, Tamil Nadu (2005); Kamrup, 
Assam (2007); Nadiad, Gujarat (2005) and Panna, Mad- 
hya Pradesh (2006). Details of epidemiological and geo- 
graphical information about study sites are reported 
elsewhere [26]. Finger prick blood was spotted on auto- 
claved Whatman filter paper strips (Number 3) from the 
symptomatic patients in active case detection surveys as 
well as from patient attending clinics. A total of 100 
microscopically-diagnosed P. vivax positive blood sam- 
ples were spotted and dried blood spots were stored at 
4.0°C. Only P. vivax infected samples were included in 
this study. This study was approved by the ethics com- 
mittee of the National Institute of Malaria Research, 
New Delhi. All blood spots were collected only after 
obtaining consent of the patients. 

DNA extraction, PGR, and DNA sequencing 

Genomic DNA was extracted from blood spots using 
QIAamp mini DNA kit (Qiagen, Germany) according to 
manufacturer instructions. Genomic DNA was eluted in 
120.0 (il triple sterile water and store in -20°C until use. 
One step modified PCR strategy was employed for 



amplification of dhfr/dhps using nested primers reported 
earlier [27,28]. The modified PCR conditions for dhfrw^s- 
initial denaturation at 95°C/5.0 minute, denaturation at 
95°C/30 second, annealing at 64°C/30 second, and exten- 
sion at 72°C/60 second for 40 cycle, and a final extension 
at 72°C/5.0 minute. The modified PCR amplification for 
dhps was same as dhfr except annealing step, which is 
55°C/30 second. PCR products were purified with Exonu- 
clease I/Shrip alkaline phosphates treatment according to 
manufacturer instruction (Fermentas, USA). Purified PCR 
products were outsourced to Macrogen Inc, Korea for 
DNA sequencing [29]. Each sample was sequenced with 
both forward and reverse primers. DNA sequences were 
edited and aligned (ClustalW method) with EditSeq and 
MegAlign module of DNA Lasergene software version 7.0 
(Madison, USA). Samples that show novel point mutations 
were re-sequenced from another independent PCR pro- 
duct. All sequences have been submitted to the GenBank 
(EU149665-EU149764, EU145878-EU145947). 

Collection of dhfr mutations data and statistical analysis 

Dhfr point mutations data from Panna, Madhya Pradesh 
(N = 30); Nadiad, Gujarat (N = 19); Mukherji Nagar, 
Delhi (N = 68); Cuttack, Orissa (N = 16); Mohan Nagar 
(N = 49) Aligarh (N = 38) and Mirzapur (N = 8), Uttar 
Pradesh; Chennai, Tamil Nadu (N = 67) Car Nicobar, 
Andaman Nicobar (N = 38); Panjim, Goa (N = 76); 
Kamrup, Assam (N = 24) and Navi Mumbai, Maharash- 
tra (N = 24) were collected from the earlier studies 
[23-25,30]. On the basis of mutations in dhfr and 
geographical location, grouping between different geo- 
graphical regions was done. Geographical regions were 
grouped according to the higher prevalence of wild type, 
double mutant, and quadruple mutant dhfr genotypes. 
Chi square test was done to check the significance of 
prevalence of dhfr mutant genotypes in geographical 
locations. Confidence interval at 95% was done to check 
the proportion dhfr genotypes in each group. 

Results 

Mutations analysis in dhfr and dhps 

For dhfr, a total of 71 P. vivax isolates were successfully 
PCR amplified and sequenced. PCR amplified dhfr was 
711 bp in length and covers complete coding region. All 
sequences were compared with wild type reference 
sequence (GenBank accession no. X98123) to detect 
point mutations in dhfr. Point mutations observed at 
dhfr were given in table 1. Approximately 50% (35/71) 
sequences at dhfr were wild type. DNA sequence align- 
ment revealed point mutations at 57(F^L), 58(S^R), 61 
(T^M), 64(V^L), 117(S^T), and 173(I^F). Point 
mutations at codon 64 (V^L) and 173 (I^F) were first 
time observed in the Indian isolates. One isolate showed 
amino acid substitution (57F^L) by changing nucleotide 
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Table 1 Distribution of Plasmodium vivax dhfr/dhps point mutations in field isolates 



DHFR DHPS 



Regions 


Sample size 


Wild type 


F57L 


S58R 


T61M 


V64L S117N/T 


I173F T69C 


C134T 


Sample size 


Wild type 


A383G 


A553G 


Nadiad 


3 


2 




1 




1 






20 


20 






Delhi 


11 


10 








1 


2 




20 


17 


3 


3 


Panna 


17 


13 




4 




3 


8 


2 


20 


20 






Chennai 


19 


3 


1 


16 


1 


16 


1 




20 


9 


11 


8 


Kamrup 


21 


7 


14 


13 


13 


13 


1 




20 


20 






Total 


71 


35 


15 


34 


14 


1 33 


1 11 


2 


100 


86 


14 


11 



Boldface: synonymous substitution. 
Values indicate number of isolates. 



TTG rather than TTA (TTC wild type). Dhfr sequence 
analysis revealed two synonymous mutations at codon 
69Y (TAT^TAC) and 134 V (GTC^GTT). Synon- 
ymous mutation 134 V (GTC^GTT) is novel and 
observed only in the wild type dhfr genotype. The 
majority of the mutated codons were S58R (47.88%), 
S117N/T (46.47%), F57L (21.12%), and T61M (19.71%). 
Study isolates revealed single (S58R), double (S58R/ 
S117N) and quadruple (F57L/S58R/T61M/S117T/) 
mutations. 

For dhps, a total of 100 P. vivax isolates were success- 
fully PGR amplified (950 bp) which covers hotspot 
region for drug resistance related point mutations 
(codon 350 to 652). All sequences were compared with 
wild type reference sequence (GenBank accession no. 
AY186730) to detect point mutations in dhps. Sequence 
analysis revealed point mutation at codon A383G and 
A583G (Table 1). Majority of the sequences were wild 
type (86%). Only single (A383G) and double (A383G/ 
A583G) point mutations were observed at dhps in the 
study isolates (Table 1). 

Geographical distribution of dhfr and dhps genotypes 

Four distinct genotypes of dhfr were observed and 
these are wild type, single mutant, double mutant and 
quadruple mutant. In total, among mutant genotypes, 
frequency of double mutant (28.16%) and quadruple 
mutant (19.71%) was higher compared with single 
mutant (2.81%). Isolates of Panna, Delhi, and Nadiad 
were dominated by wild type genotype (Table 2). 
Frequency of mutant dhfr genotypes was highest in 
isolates of Chennai (84.24%) and Kamrup (66.67%) 
regions. In contrast, isolates from Delhi (9.0%), Panna 
(23.52%), and Nadiad (33.34%) showed relatively low 
proportion of mutant genotypes. Geographical distribu- 
tion of various genotypes observed at dhfr is given in 
Table 2. Double mutant genotype was observed to be 
predominated in Chennai isolates (78.94%) whereas 
quadruple mutant genotype was dominated in Kamrup 
region (61.90%). The silent point mutation T69C was 
observed in wild type dhfr sequences only. The point 



mutation T69C was found at a high frequency in iso- 
lates of Panna (47.05%) but in lesser extent in Delhi 
(18.18%) and Chennai (5.26%) isolates. Another silent 
point mutation C134T was limited to isolates of Panna 
(11.76%). 

Three distinct dhps genotypes were observed in the 
field isolates: wild type, single mutant and double 
mutant. The wild type genotype was observed in high 
frequency in field isolates of Delhi (85%) and Chennai 
(55%) whereas isolates from Panna, Nadiad, and Kamrup 
were exclusively wild type. The mutant genotypes (single 
and double) were only observed in isolates from Delhi 
and Chennai (Table 2). 

Tandem repeat variation, distribution, and association 
with mutations 

At dhfry four tandem repeat variants were observed and 
designated as Type 1-4, on the basis of six amino acid 
deletion/insertion (Figure 1). In total, tandem repeat 
Type 2 (84.50%) was the major repeat polymorph 
observed in field isolates (Table 3). Type 2 tandem 
repeat was observed in high frequency in all study sites 
except for Delhi where Type 1 predominates. Distribu- 
tion of tandem repeat variants in different study sites is 
given in Table 3. Tandem repeat variant was mono- 
morphic (Type-2) in field isolates of Nadiad, Chennai 
and Kamrup, whereas it was polymorphic in field iso- 
lates of Delhi and Panna. This suggests that 

Table 2 Region wise distribution of Plasmodium vivax 
dhfr and dhps genotypes 



DHFR (%) DHPS (%) 



Region 


Sample 
size 


WT 


SM 


DM 


QM 


Sample 
size 


WT 


SM 


DM 


Nadiad 


3 


66.67 




33.33 




20 


100 






Delhi 


11 


90.9 


9.1 






20 


85 




15 


Panna 


17 


76.47 


5.88 


17.65 




20 


100 






Chennai 


19 


15.79 




78.95 


5.26 


20 


45 


15 


40 


Kamrup 


21 


33.34 




4.76 


61.9 


20 


100 






Total 


71 


49.29 


2.81 


28.16 


19.71 


100 


86 


3 


11 



Prajapati et al. Malaria Journal 201 1, 10:102 
http://www.malariajournal.eom/content/10/1/102 



Page 4 of 7 



256 QGGGDNTSGGDNTHGGDNTHGGDNADKLQT Type 1 

256 QGGGDNTSGGDNTHG GDNADKLQT Type 2 

256 QGGGDNT HGGDNADKLQT Type 3 

256 QGGGDNT SGGDNADKLQT Type 4 

Figure 1 Tandem repeat variation in Plasmodium vivax DHFR. 



f ^ 

754 GEAKLTNGEGKLTNGEAKLTNGEGKLTNGEAKLTNGEGKLTNGDAKLTNGDSKLTNG Type A 

754 GEAKLTNGEGKLTNGEAKLTNGEGKLTNG DSKLTNG Type B 

754 GEAKLTNGEGKLTNGEAKLTNGEGKLTN GDAKLTNGDSKLTNG Type C 

754 GEAKLTNGEGKLTNGEAKLTNGEGKLTNGEAKLTNGEGKLTN GDSKLTNG Type D 

754 GEAKLTNGEGKLTNG DSKLTNG Type E 

754 GEAKLTNGEGKLTNGEAKLTN GEAKLTNGEGKLTNGDAKLTNGDSKLTNG Type F 

754 GEAKLTNGEGKLTNGE GKLTNGEAKLTNGEGKLTN GDSKLTNG Type G 

754 GEAKLTNGEGKLTNGEGKLTNGEGKLTNGEAKLTNG DSKLTNG Type H 

Figure 2 Tandem repeat variation in Plasmodium vivax DHPS. 



monomorphic nature of tandem repeat was only 
observed in a region where mutant dhfr genotypes pre- 
dominate. In study isolates, all mutant dhfr genotypes 
were observed in isolates having Tandem repeat Type-2. 

At dhps, eight tandem repeats variants were observed 
and designated as Type A-H on the basis of deletion/ 
insertion of seven amino acid repeat whereas (Figure 2). 
Tandem repeat Type-B was the major variant in the 
field isolates (75.25%) and predominate in all study sites 
(Table 3). Distribution of tandem repeat variants in dif- 
ferent study sites is given in Table 3. Seven tandem 
repeat variants were found in field isolates of Kamrup, 
three in Panna, Chennai, and Nadiad and only two in 
Delhi. Point mutations were only observed in field iso- 
lates having tandem repeat Type B. 

Progressive increase in frequency of dhfr quadruple 
mutant 

In Kamrup region, dhfr point mutation data from 2002- 
2003 revealed higher frequency of double mutant 
(54.16%) and low prevalence of wild (29.16%) and quad- 
ruple mutant (16.67%) genotypes [23]. Present study 



Table 3 Distribution of tandem repeat variants at DHFR 
and DHPS in field isolates 



DHFR 






Regions (%) 




Total N 

(%) 


Repeat 
Type 


Deliii 


Nadiad Panna 


Chiennai 


Kamrup 




Type 1 


54.54 


0 


11.76 


0 


0 


8 (11.26) 


Type 2 


36.34 


100 


76.47 


100 


100 


60 (84.5) 


Type 3 


9.09 


0 


0 


0 


0 


1 (1.4) 


Type 4 


0 


0 


11.76 


0 


0 


2 (2.81) 


Total (N) 


11 


3 


17 


19 


21 


71 


DHPS 














Type A 


10 


10 


0 


5 


21.05 


9.27 


Type B 


90 


85 


83.34 


90 


26.31 


75.25 


Type C 


0 


0 


5.55 


5 


0 




Type D 


0 


5 


5.55 


0 


15.79 


5.15 


Type E 


0 


0 


0 


0 


10.52 


2.06 


Type F 


0 


0 


0 


0 


10.52 


2.06 


Type G 


0 


0 


5.55 


0 


10.52 


3.09 


Type H 


0 


0 


0 


0 


5.26 


1.03 


Total (N) 


20 


20 


18* 


20 


19* 


97 



*: DNA sequences of tandem repeat region was not good in some isolates 
and therefore not included in analysis of tandem repeats. 



conducted in 2006-2007 revealed quite surprising result 
with respect to the data of dhfr point mutations studied 
in 2002-03. In the present study a high prevalence of 
the quadruple mutant genotype (61.90%) and an almost 
absence of the double mutant genotype (4.16%) were 
observed. However, the lower prevalence of wild type 
remains the same. The interesting finding was the dras- 
tic increase in the frequency of quadruple mutant dhfr 
genotype in Kamrup region (p = <0.0001) within a five 
year (2002 to 2007) span. This suggests that there is 
progressive increase in the frequency of quadruple 
mutant dhfr genotype. 

Geographical clustering of dhfr genotypes among Indian 
continent 

Dhfr point mutations data of earlier published work (n = 
477) revealed single to quadruple mutant genotypes with 
their varied prevalence in different geographical regions 
of India (Figure 3). Comprehensive analysis of dhfr 
point mutations revealed three distinct pattern which 
are- 1) high prevalence of wild type, 2) high prevalence 
of double mutant type, and 3) high prevalence of quad- 
ruple mutant type. On the basis of above categorization, 
study sites were grouped together to see any trend of 
significance between geographical location and preva- 
lence of dhfr genotypes. Three distinct geographical 
clusters were appeared at significant level (p < 0.0001). 
These geographical clusters are- Type-1) covers geogra- 
phical areas of North India and majority of isolates were 
wild type (81.64%, 95% CI = ± 4.9) and very low propor- 
tion of single and or double mutant. Type-2) covers 
geographical areas of South India and majority of iso- 
lates were double mutant type (80.54%, 95% CI = ± 
6.22) and very low proportion of wild and single mutant 
types. Type-3) covers geographical areas of north-east- 
ern state and Island region where quadruple mutant 
(42.68%, 95% CI = ± 10.72) predominate over wild 
(29.26%) and double mutant type (28.04%). Details of 
geographical regions falling in specific cluster are given 
in figure 3. Comparison of Type-1 cluster with Type-2 
or mutant versus wild types (p < 0.0001), and Type-2 
cluster with Type-3 or double mutant versus quadruple 
mutant (p < 0.0001) strongly supports grouping or clus- 
tering of wild, double mutant, and quadruple mutant 
dhfr genotypes in three distinct geographical regions of 
the country. 
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1: Delhi, 2: Mohan Nagar, 3: Aligarh 
4: Mirzapur, 5: Panna, 6: Nadiad 
7: Navi Mumbai, 8: Goa, 9: Chennai 
JO: Cuttack, 11: Kamrup, 12: Car Nicobar 



.11 




O Type 3 



Wild type 
Single mutanlT" 
Double mutanifi| 

I Triple mutant 

I Quadruple mutant 



Figure 3 Frequencies of Plasmodium vivax dhfr genotypes and their geographical clustering among Indian subcontinent. 



Discussion 

Due to emergence and widespread of drug resistant 
strains of human Plasmodium species against che- 
motherapeutic intervention agents, clinical trial for 
monitoring drug efficacy over regular intervals as well as 
monitoring mutation data in relevant genes related with 
anti-malarial resistance are the crucial steps in deciding 
scope of a drug in field. This study uncovers 1) progres- 
sive increase in the frequency of quadruple mutant dhfr 
genotypes, 2) limited point mutations in dhps among 



Indian isolates, and 3) distinct pattern of anti-folate 
resistance in Indian sub-continent linked with geogra- 
phical regions. 

Plasmodium vivax is still susceptible to chloroquine in 
India [30,31] and anti-folate drugs are not used for treat- 
ment of vivax malaria, therefore, the appearance of point 
mutations in dhfr/dhps is surprising. This has been 
explained earlier because of the use of sulphadoxine-pyri- 
methamine to treat chloroquine-resistant P, falciparum is 
creating selection pressure in the P. vivax population. 
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This has been reflected in the isolates of areas with 
sympatricity of P, falciparum and P, vivax [23,25], which 
presented a higher proportion of double and quadruple 
mutant dhfr genotypes than other regions. Cotrimoxazole 
(sulfamethoxazole and trimethoprim), an anti-microbial 
drug that target bacterial dhfr/dhps, is widely used in 
bacterial infection as prophylaxis, and has potential to 
select mutant dhfr/dhps of malaria parasites [32,33]. 
Cotrimoxazole is being used in India to treat bacterial 
infection in wide spectrum [34,35], which suggests that it 
could be a potential factor to select anti-folate point 
mutations in the Indian subcontinent. Other factors, i.e. 
presumptive treatment of malaria without diagnosis and 
use of S/P by private practitioners, cannot be ruled out. 

The comprehensive analysis of point mutations in dhfr 
among Indian sub-continent revealed three distinct clus- 
ters/groups/types linked with geographical locations. 
North and south India revealed a prevalence of wild and 
double mutant dhfr^ respectively, whereas north-eastern 
and island regions revealed a high prevalence of quadru- 
ple mutant dhfr genotype. Single to quintuple mutations 
have been reported in dhfr but only triple to quintuple 
mutant alleles are known to confer high level of resis- 
tance [36-38]. This strongly suggests that P, vivax iso- 
lates from northern and southern part of India could be 
sensitive to anti-folate drug, whereas isolates from 
north-eastern state and Island regions will be resistant 
to anti-folate drug. 

A significant association between dhfr mutations and 
tandem repeat polymorphism was observed in studied iso- 
lates. The mutated dhfr genotype was exclusively asso- 
ciated with Type-2 tandem repeat and is strongly 
supported by previous study from India [23]. Further 
quadruple mutant dhfr alleles were exclusively associated 
with type 1 tandem repeat and are supported by earlier 
studies from India [23], Thailand [38] and Myanmar [39], 
however exception to this association was only showed by 
a single isolate from Myanmar [39]. These observations 
speculate dhfr allele with type 1 tandem repeat could be 
more prone for mutations and development of higher 
level of resistance conferring genotypes. Therefore, it 
could act as molecular marker to predict the risk of 
mutant/higher level resistance conferring mutant geno- 
types in any geographical area of Indian subcontinent. 

Frequency of quadruple mutant dhfr genotype was 
progressively increased in Kamrup isolates from 2005 to 
2007 [23]. The high proportion of quadruple mutant 
dhfr genotype in endemic area may be: 1) because of the 
characteristics of mutant parasites associated with 
higher gametocytogenesis [40-43] which favour trans- 
mission of mutant parasites compared to wild type dhfr, 
and 2) north-eastern states of India are highly endemic 
for malaria and international borders with neighbouring 
countries surround the region. The borders are very 



porous and illegal migrations of people across the bor- 
ders are very common thus, importing drug resistant 
strains. The higher frequency of quadruple mutant dhfr 
genotypes in the isolates of Myanmar (71%) and Thai- 
land (96%) [39,44], supports for inflow of quadruple 
mutant dhfr genotype in north-eastern states of India. 

The development and spread of drug resistant parasite 
strains is a major obstacle to the malaria elimination pro- 
gramme. Therefore, it is essential to identify drug resis- 
tance areas/regions on the basis of point mutations in 
order to manage the national anti-malarial drug policy. In 
this regard, this study revealed a higher and increasing 
prevalence of quadruple mutant dhfr genotype in north- 
eastern and Island regions, therefore, additional caution 
may be taken for treatment of vivax malaria in these 
regions to stop the flourishing of quadruple mutant dhfr 
genotypes in remaining part of country. The study con- 
cluded that P, vivax isolates from northern and southern 
part of India could be susceptible to SP except for the iso- 
lates from Car Nicobar (Island) and Assam (north-eastern 
region). Geographical clustering of dhfr mutant genotypes 
would provide a rationale for a more appropriate national 
anti-malarial drug policy. 



Acknowledgements 

This work is supported by Department of Biotechnology, New Delhi. SKP is 
an ICMR Postdoctoral Fellow. The authors are thankful to the patients for 
their co-operation during the study. The authors are also thankful to 
scientist/staff of National Institute of Malaria Research, New Delhi as well as 
from field station for their support and help during the study. Authors 
expressed their deep gratitude to Dr. Hema Joshi, who passed away recently 
and she had contributed significantly to the present study. 

Author details 

^Genetics and Molecular Biology Laboratory, National Institute of Malaria 
Research (NIMR), Sector-8, Dwarka, New Delhi 110077, India. ^NIMR field 
station, Assam, India. 

Authors' contributions 

SKP; Experiment design, experimental work, data analysis, manuscript writing 
HJ: Conceptual design of work, data analysis, manuscript writing, 
VD: Sample collection, parasite identification, and experimental design 
VKD; Conceptual design of work. Overall supervision of work, and 
manuscript writing 

All authors read and approved the final manuscript. 
Competing interests 

The authors declare that they have no competing interests. 

Received: 13 January 201 1 Accepted: 24 April 201 1 
Published: 24 April 201 1 

References 

1. WHO: Annual Report. World Health Organization Geneva; 2008. 

2. Andrade BB, Reis-Filho A, Souza-Neto SM, Clarencio J, Camargo LM, Barral A, 
Barral-Netto M: Severe Plasmodium vivax malaria exhibits marked 
inflammatory imbalance. Malar J 2010, 9:13. 

3. Kochar DK, Das A, Kochar SK, Saxena V, Sirohi P, Garg S, Kochar A, Khatri MP, 
Gupta V: Severe Plasmodium vivax malaria: a report on serial cases from 
Bikaner in northwestern India. Ann J Trop Med Hyg 2009, 80:194-198. 

4. Kochar DK, Saxena V, Singh N, Kochar SK, Kumar SV, Das A: Plasmodium 
vivax malaria. Emerg Infect Dis 2005, 11:132-134 



Prajapati et al. Malaria Journal 201 1, 10:102 
http://www.malariajournal.eom/content/10/1/102 



5. Genton B, D'Acremont V, Rare L, Baea K, Reeder JC, Alpers MP, Muller I: 
Plasmodium vivax and mixed infections are associated with severe 
malaria in children: a prospective cohort study from Papua New Guinea. 
PLoS Med 2008, 5:el27. 

6. Rogerson SJ, Carter R: Severe vivax malaria: newly recognised or 
rediscovered. PLoS Med 2008, 5:el36. 

7. Tjitra E, Anstey NM, Sugiarto P, Warikar N, Kenangalem E, Karyana M, 
Lampah DA, Price RN: Multidrug-resistant Plasmodium vivax associated 
with severe and fatal malaria: a prospective study in Papua, Indonesia. 
PLoS Med 2008, 5:el28. 

8. Chulay JD, Watkins WM, Sixsmith DG: Synergistic antimalarial activity of 
pyrimethamine and sulfadoxine against Plasmodium falciparum in vitro. 
Am J Trop Med Hyg 1984, 33:325-330. 

9. Cowman AF, Morry MJ, Biggs BA, Cross GA, Foote SJ: Amino acid changes 
linked to pyrimethamine resistance in the dihydrofolate reductase- 
thymidylate synthase gene of Plasmodium falciparum. Proc Natl Acad Sci 
USA 1988, 85:9109-9113. 

10. Triglia T, Meriting JG, Wilson C, Cowman AF: Mutations in dihydropteroate 
synthase are responsible for sulfone and sulfonamide resistance in 
Plasmodium falciparum. Proc Natl Acad Sci USA 1997, 94:13944-13949. 

11. Sehgal PN, Sharma MID, I SS, Gogal S: Resistance to chloroquine in 
falciparum malaria in Assam state, India. J Commun Dis 1973, 5:175-180. 

12. National Vector Borne Disease Control Programme. [http://www.nvbdcp. 
gov.in]. 

13. Rieckmann KH, Davis DR, Hutton DC: Plasmodium vivax resistance to 
chloroquine? Lancet 1989, 2:1183-1184. 

14. Garavelli PL, Corti E: Chloroquine resistance in Plasmodium vivax: the first 
case in Brazil. Trans R Soc Trop Med Liyg 1992, 86:128. 

15. Garg M, Gopinathan N, Bodhe P, Kshirsagar NA: Vivax malaria resistant to 
chloroquine: case reports from Bombay. Trans R Soc Trop Med Liyg 1995, 
89:656-657. 

16. Murphy GS, Basri H, Purnomo , Andersen EM, Bangs MJ, Mount DL, Gorden J, 
Lai AA, Purwokusumo AR, Harjosuwarno S, et ah Vivax malaria resistant to 
treatment and prophylaxis with chloroquine. Lancet 1993, 341:96-100. 

17. Myat Phone K, Myint 0, Myint L, Thaw Z, Kyin HIa A, Nwe Nwe Y: 
Emergence of chloroquine-resistant Plasmodium vivax in Myanmar 
(Burma). Trans R Soc Trop Med Liyg 1 993, 87:687. 

18. Schuurkamp GJ, Spicer PE, Kereu RK, Bulungol PK, Rieckmann KH: 
Chloroquine-resistant Plasmodium vivax in Papua New Guinea. Trans R 
Soc Trop Med Hyg 1992, 86:121-122. 

19. Schwartz IK, Lackritz EM, Patchen LC: Chloroquine-resistant Plasmodium 
vivax from Indonesia. N Engl J Med 1991, 324:927. 

20. Hasugian AR, Tjitra E, Ratcliff A, Siswantoro H, Kenangalem E, Wuwung RM, 
Purba HL, Piera KA, Chalfien F, Marfurt J, Penttinen PM, Russell B, 

Anstey NM, Price RN: In vivo and in vitro efficacy of amodiaquine 
monotherapy for treatment of infection by chloroquine-resistant 
Plasmodium vivax. Antlmlcrob Agents Ctiennottier 2009, 53:1094-1099. 

21. Russell B, Chalfein F, Prasetyorini B, Kenangalem E, Piera K, Suwanarusk R, 
Brockman A, Prayoga P, Sugiarto P, Cheng Q, Tjitra E, Anstey NM, Price RN: 
Determinants of in vitro drug susceptibility testing of Plasmodium vivax. 
Antlmlcrob Agents CtiemotLier 2008, 52:1040-1045. 

22. Kumar A, Valecha N, Jain T, Dash AP: Burden of malaria in India: 
retrospective and prospective view. Am J Trop Med Liyg 2007, 77:69-78. 

23. Alam MT, Bora H, Bharti PK, Saifi MA, Das MK, Dev V, Kumar A, Singh N, 
Dash AP, Das B, Wajihullah , Sharma YD: Similar trends of pyrimethamine 
resistance-associated mutations in Plasmodium vivax and P. falciparum. 
Antlmlcrob Agents ChemotLier 2007, 51:857-863. 

24. Kaur S, Prajapati SK, Kalyanaraman K, Mohmmed A, Joshi H, Chauhan VS: 
Plasmodium vivax dihydrofolate reductase point mutations from the 
Indian subcontinent. Acta Trop 2006, 97:174-180. 

25. Prajapati SK, Joshi H, Valecha N, Reetha AM, Eapen A, Kumar A, Das MK, 
Yadav RS, Rizvi MA, Dash AP: Allelic polymorphism in the Plasmodium 
vivax dihydrofolate reductase gene among Indian field isolates. Clin 
Microbiol Infect 2007, 13:331-334. 

26. Prajapati SK, Verma A, Adak T, Yadav RS, Kumar A, Eapen A, Das MK, 
Singh N, Sharma SK, Rizvi MA, Dash AP, Joshi H: Allelic dimorphism of 
Plasmodium vivax gam-1 in the Indian subcontinent. Malar J 2006, 5:90. 

27. Imwong M, Pukrittakayamee S, Looareesuwan S, Pasvol G, Poirreiz J, 
White NJ, Snounou G: Association of genetic mutations in Plasmodium 
vivax dhfr with resistance to sulfadoxine-pyrimethamine: geographical 
and clinical correlates. Antlmlcrob Agents CLiemotLier 2001, 45:3122-3127. 



Page 7 of 7 



28. Menegon M, Majori G, Severini C: Genetic variations of the Plasmodium 
vivax dihydropteroate synthase gene. Acta Trop 2006, 98:196-199. 

29. Macrogen Inc: Korea.[http://www.macrogen.com]. 

30. Valecha N, Joshi H, Eapen A, Ravinderan J, Kumar A, Prajapati SK, Ringwald P: 
Therapeutic efficacy of chloroquine in Plasmodium vivax from areas with 
different epidemiological patterns in India and their Pvdhfr gene 
mutation pattern. Trans R Soc Trop Med Liyg 2006, 100:831-837. 

31. Nandy A, Addy M, Maji AK, Bandyopadhyay AK: Monitoring the 
chloroquine sensitivity of Plasmodium vivax from Calcutta and Orissa, 
India. Ann Trop Med Parasltol 2003, 97:215-220. 

32. Jelinek T, Kilian AH, Curtis J, Duraisingh MT, Kabagambe G, von 
Sonnenburg F, Warhurst DC: Plasmodium falciparum: selection of serine 
108 of dihydrofolate reductase during treatment of uncomplicated 
malaria with co-trimoxazole in Ugandan children. Am J Trop Med Liyg 
1999, 61:125-130. 

33. Kiara SM, Okombo J, Masseno V, Mwai L, Ochola I, Borrmann S, Nzila A: In 
vitro activity of antifolate and polymorphism in dihydrofolate reductase 
of Plasmodium falciparum isolates from the Kenyan coast: emergence of 
parasites with lle-164-Leu mutation. Antlmlcrob Agents CLiemotLier 2009, 
53:3793-3798. 

34. Raizada N, Chauhan LS, Babu BS, Thakur R, Khera A, Wares DF, Sahu S, 
Bachani D, Rewari BB, Dewan PK: Linking HIV-infected TB patients to 
cotrimoxazole prophylaxis and antiretroviral treatment in India. PLoS One 

2009, 4:e5999. 

35. Bansa V, Jose V, Medhi B, Pandhi P: Decrease Use of Cotrimoxazole in 
Clinical Practice: Is it Justified? JK SCIENCE 2010, 12:214-215. 

36. Hastings MD, Maguire JD, Bangs MJ, Zimmerman PA, Reeder JC, Baird JK, 
Sibley CH: Novel Plasmodium vivax dhfr alleles from the Indonesian 
Archipelago and Papua New Guinea: association with pyrimethamine 
resistance determined by a Saccharomyces cerevisiae expression 
system. Antlmlcrob Agents CLiemotLier 2005, 49:733-740. 

37. Hawkins VN, Joshi H, Rungsihirunrat K, Na-Bangchang K, Sibley CH: 
Antifolates can have a role in the treatment of Plasmodium vivax. Trends 
Parasltol 2007, 23:213-222. 

38. Imwong M, Pukrittayakamee S, Renia L, Letourneur F, Charlieu JP, 
Leartsakulpanich U, Looareesuwan S, White NJ, Snounou G: Novel point 
mutations in the dihydrofolate reductase gene of Plasmodium vivax: 
evidence for sequential selection by drug pressure. Antlmlcrob Agents 
CLiemotLier 2003, 47:1514-1521. 

39. Na BK, Lee HW, Moon SU, In TS, Lin K, Maung M, Chung GT, Lee JK, Kim TS, 
Kong Y: Genetic variations of the dihydrofolate reductase gene of 
Plasmodium vivax in Mandalay Division, Myanmar. Parasltol Res 2005, 
96:321-325. 

40. Mockenhaupt FP, Teun Bousema J, Eggelte TA, Schreiber J, Ehrhardt S, 
Wassilew N, Otchwemah RN, Sauerwein RW, Bienzle U: Plasmodium 
falciparum dhfr but not dhps mutations associated with sulphadoxine- 
pyrimethamine treatment failure and gametocyte carriage in northern 
Ghana. Trop Med Int HealtLi 2005, 10:901-908. 

41. Adedeji AA, Fehintola FA, Fateye BA, Happi TC, Amoo AO, Gbotosho GO, 
Sowunmi A: Plasmodium falciparum malaria in Nigerian children during 
high and low transmission seasons: gametocyte carriage and response 
to oral chloroquine. J Trop Pedlatr 2005, 51:288-294. 

42. Sowunmi A, Fateye BA, Adedeji AA, Fehintola FA, Gbotosho GO, Happi TC, 
Oduola AM: Predictors of the failure of treatment with chloroquine plus 
chlorpheniramine, in children with acute, uncomplicated, Plasmodium 
falciparum malaria. Ann Trop Med Parasltol 2005, 99:331-338. 

43. Tjitra E, Suprianto S, Anstey NM: Higher gametocyte prevalence following 
failure of treatment of Plasmodium falciparum malaria with sulfadoxine- 
pyrimethamine and the combination of chloroquine plus sulfadoxine- 
pyrimethamine: implications for progression of anti-folate resistance. 
Trans R Soc Trop Med Hyg 2002, 96:434-437. 

44. Lu F, Lim CS, Nam DH, Kim K, Lin K, Kim TS, Lee HW, Chen JH, Wang Y, 
Sattabongkot J, Han ET: Mutations in the antifolate-resistance-associated 
genes dihydrofolate reductase and dihydropteroate synthase in 
Plasmodium vivax isolates from malaria-endemic countries. Am J Trop 
Med Hyg im, 83:474-479. 



doi:1 0.1 1 86/1 475-2875-1 0-1 02 

Cite this article as: Prajapati et al:. Molecular epidemiology of 
Plasmodium vivax anti-folate resistance in India. Malaria Journal 201 1 
10:102. 



